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We demonstrate that Majorana neutrinos can form Cooper pairs due to long-range attractive 
forces and show BCS superfluidity in a class of mass varying neutrino dark energy models. We 
describe the condensates for Majorana neutrinos and estimate the value of the gap, critical temper- 
ature and Pippard coherence length for a simple neutrino dark energy model. In the strong coupling 
regime bosonic degree of freedom can become important and Bose- Einstein condensate may govern 
the dynamics for the mass varying neutrino models. Formation of the condensates can significantly 
alter the instability scenario in the mass varying neutrino models. 
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Some time back neutrino superfluidity was studied for 
Dirac neutrinos [lj], in which the left-handed and the 
right-handed neutrinos form Cooper pairs due to the at- 
tractive force originating from their Yukawa interactions 
with the Higgs scalar. Unfortunately this interesting con- 
cept could not be applied to any realistic situation in 
astrophysics or cosmology. We extend this formalism 
for Majorana neutrinos and show that the superfluid- 
ity of relic neutrinos could be important when one con- 
siders their interactions with very light scalar field, like 
quintessence. We study the superfluidity of Majorana 
neutrinos in the context of the mass varying neutrinos 
(MaVaNs) models JIB and show that it can solve the 
stability problem [J, [5[ of the MaVaNs naturally. 

The interaction between the neutrinos in a MaVaN sce- 
nario is known to be attractively, [7|, la] due to the pres- 
ence of a quintessence field called acceleron. It is well 
known that these models are unstable when the neutri- 
nos become non-relativistic i.e. their pressure p v « 
0, Sj- The instability saturates when the degeneracy 
pressure balances the attractive force and the final state 
can evolve as ACDM. It should be noted that this insta- 
bility does not arise in a certain class of models involving 
super-acceleration 0] 

We show that this attractive interaction can lead to 
neutrino superfluidity in MaVaNs by formation of the 
Cooper pairs. If the size of the Cooper pairs is smaller 
than the length scales relevent for the dark energy dy- 
namics, the dynamics of a scalar field describing the Bosc- 
Einstein condensate could be applied for studying the 
evolution of the system. The inclusion of the conden- 
sate dynamics alters the instability scenario significantly. 
Firstly, there would be no degeneracy pressure in the 
bosonic system. Moreover the coupling between the neu- 
trinos and the scalar field will be changed. However the 
attractive force, if any, between the condensates and the 
acceleron should be balanced by the Heisenberg uncer- 
tainty. This kind of stable structures are known in the 
literature as Boson stars . Next we demonstrate that 
the stability calculations considered earlier 0, [B| are al- 
tered and the new stability criteria can be satisfied by 
the different models of dark energy potentials. Conden- 
sates with Majorana neutrinos has been discussed in the 



literature [13], but here we develope a statistical formal- 
ism following ref. [lj. 

The first problem one encounters while dealing with 
the Majorana neutrino supcfluidity is the chemical poten- 
tial. Since the Majorana particles self annihilate, num- 
ber of particles is not conserved, and hence, their number 
operator and the chemical potential vanishes in equilib- 
rium. However, in the early universe, these are not of 
any concern: the Majorana neutrinos have two helicity 
components which can vary with time according to the 
helicity-flip rate. One can define the chemical potential 
to the extent the helicity is conserved. Typically the ratio 
of the helicity-flip rate to the current Hubble expansion 
H rates is 



G 2 T 3 gr 



10" 



(1) 



where g is the effective degrees of freedom relevent for the 
temperature T [l]| . Thus the helicity- flipping rate of the 
Majorana particles ceases below a certain temperature, 
when the particles move apart from each other due to 
the expansion of the universe at a faster rate compared 
to their self annihilation rate [3]. 

Another important problem is to show how the freely 
streaming streaming neutrinos can become degenerate 
and exhibit superfluidity: In any MaVaN model there 
is an interaction between the scalar field and the neutri- 
nos. Dynamics of the neutrinos can be described by the 
following kinetic equation [Afshordi et.al. in Ref. 



drj 



07 1 Vm t 



dp 



= 



(2) 



where a and rj are the scale factor and the conformal time, 
defined in Robertson- Walker metric, respectively. The 
last term includes the effective neutrino mass variation 
due to the scalar field. In the absence of the last term 
on the left hand side, / is given by the usual Fermi-Dirac 
distribution. However, when small perturbations of the 
type Sf = A(p)exp[i(k ■ x — wrj\ and 8m v = T,exp[i(k ■ 
x — tor]] are considered, the system becomes unstable and 
eq.(l) gives, in a sub-Hubble regime, 



wA(p) = k • u — 7 



dp 
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Only Lo/k — c s appear in the equation and there is no 
preferred scale in this equation. c s can be found (Af- 
shordi et. al. in Ref.0]) to be ±y/—l in a nonrelativisitic 
regime. Thus the instability that grows at smaller length 
scale. Since there is no preferred length scale in equation 
([3]), the instability can continue to grow on the smaller 
length scales. However in a realistic situation the in- 
stability may saturate when the deg eneracy pressure by 
the neutrinos become important 18]. At this stage the 
attractive interaction between the degenerate neutrinos 
may induce the phenomenon of superfluidity. 

Consider a Majorana mass term of the left-handed neu- 
trinos in a model that provides an attractive long-range 
force due to the exchange of a scalar field <j>\ 



Cm = v M K + U 0] u m 



(4) 



Although we present the formalism for one gneration, it 
can be easily extended to the realistic case with three 
neutrinos. 

The Majorana field vm, defined in terms of the left- 
handed neutrinos vl and its CP conjugate field v c r: 



vm = v L + \v c R 
satisfies the condition: 

V M = K + ^ C rT = A*fM- 



(5) 
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where A is the Majorana phase, |A| 2 = 1. We work in the 
Weyl representation: 
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with a* = [I 2 ,<Ti\; cr M = [h,-ai] 



where 1% is a 2 x 2 unit matrix and a 1 are the Pauli ma- 
trices. In this basis, 75 = Z7°7 1 7 2 7 3 = diag(— 1%, 1%) 
is diagonal, and the left and right-handed fields become: 
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The Majorana neutrinos can now be expressed as: 



vm = 
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so that the Lagrangian density associated with the mass 
term (equation 2]) becomes 



In terms of the component fields ip this becomes 
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One of the key ingredient in theory of superconduc- 
tivity is to have an overall attractive interaction be- 
tween its particles. In the case of a metal the Coulomb 
interaction between the electrons, in Fourier space is 
Vcoui = e 2 /(k 2 +Kp), where i0D _1 is the typical shield- 
ing distance. Vcoul is always repulsive i.e. V cou i > 0. 
However the superconductivity arises as the the electrons 
in the metal also have an attractive interaction V p h aris- 
ing due to their interaction with the phonons. One can 
show that for a superconductor V p h + Vcoui < 0. Un- 
der this condition it is energetically more favourable for 
particles to form pairs. Momenta of the particles in the 
pair states are directed opposite with each other and have 
spin in the opposite directions. When the energy mini- 
mization is carried out for the wave function containing 
the pair states either occupied by the two particles or 
none the gap condition naturally arises (lBj . In the case 
of MaVaN scenario, as implied by Eqs.f MTU]) there is an 
attractive interaction between the neutrons. In fact the 
famous instability in this scenario arises In fact the in- 
stability in MaVaN scenario arises precisely when this 
interaction dominates over the gravity jij. Since there 
is no other interaction that can make up for the overall 
repulsive interaction, the above condition of supercon- 
ductivity can be satisfied. 

The gauge boson exchange would give repulsive force 
between two left-handed fields, so the only possible con- 
densate would correspond to a spin-0 pairing of the left- 
handed neutrinos with the right-handed antineutrinos: 



(ip a tp b ) = e ab D . 
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The mean field approximation would then give us the 
interaction Hamiltonian with the condensate D: 

H™ F = -2 C [A* 2 0l fa D + A 2 0-1 ^ £,♦] £ab . (12) 

We shall now express the Majorana field in terms of the 
creation and the annihilation operators as 



Cm = K + U4>]vm v m = K + 



A>V + tyH] ■ Mx) = E (fpsUpse-** + X*f^v ps e^) 



(13) 



It ought to be mentioned that the present mass of scalar 
field differ widely from 10~ 4 eV as in Fardon et. al. in 
Ref. to ?Ti0 > H e.g Bjaelde et. al. in Ref. [H]. In 
what follows, we deviate from the initial model of Fardon 
et. al. and consider the case with > H. For small 
energy and momentum transfers, interaction term can be 
written as 



The component fields are then related to the creation and 
annihilation operators through the relation 
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The interaction Hamiltonian can then be written in terms 
of the creation and annihilation operators as: 
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where e = \/p 2 + m 2 . 

In models of mass varying neutrinos the number den- 
sity of the Majorana neutrinos becomes proportional to 
the inverse of neutrino mass. In addition, the effective 
interaction Hamiltonian also does not conserve particle 
number. Since the treatment is based on grand canoni- 
cal ensemble, this requires a self-consistent treatment to 
determine when the condensates become nonvanishing. 
The complete Hamiltonian (H) is obtained by adding the 
interaction part H± and the free-particle Hamiltonian 



(15) 



One can also write the complete Hamiltonian in a so 
called standard or canonical form in which it resembles 
with the free particle Hamiltonian in Eq.(15) [HJ as 



H 



pT^pT 



(16) 



where E 2 = (e — fi) 2 + k 2 and \x is the chemical potential. 

For the late universe when the neutrino become non- 
relativistic, T <§C m, one can write its chemical potential 
following Ref.pj} as 

jti(i) = m + (/in - m)T{t)/T D 

where Tjy is the decoupling temperature and T(t) can be 
regarded as the current temperature. In the late universe 
the second term on the right hand side can be negligible 
compared to the first term. 

One can have a time-dependent transformation that 
can relate the complete Hamiltonian H = Hn + Hf IF 
with the standard form given by Eq.(16) [EEg. A rela- 
tion between the annihilation and the creation operators 
in both the Hamiltonians is given by 

b p] = cos6»e l(a+et) /pT -sin0e i(Q+et) /I pi 

b pl = cos0e^ +et Vpi+sm#e^ + ^/i pT , (17) 

with D\* 2 = \D\e 2la , tan 26 = n/(e - (i) and K = 
2 C\D\m v /e. A consistent solution for the nonvanishing 
condensate D ^= requires a = 7r/2, which has contri- 
butions from both the condensate D as well as from the 
Majorana phase A*. The magnitude of the gap is deter- 
mined by the consistency condition that the value of the 
condensate is same as that of the value obtained by the 



canonical transformation. In other words, if we express 
the condensate in terms of the density matrix (p) as 



(i>a ^b) 



(18) 



the density matrix (p) satisfies the consistency condition 
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This condition translates into 
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(20) 



whose solution gives us the magnitude of the gap. This 
integral is divergent and it should be cut off with the 
upper limit A. The condensates form due to the attrac- 
tive force between the neutrino and the scalar field in 
the MaVaN scenario. For the early times when the neu- 
trinos were in thermal contact with the other species in 
the universe, this attractive interaction may not be very 
important. Thus the values of A can be estimated from 
the energy scales below which the attractive interaction 
can be felt by the neutrinos become important. 
Solving this equation we obtain the gap 



1/3 



(21) 



where x = 2tt 2 /[CmKSTr 2 ^) 1 ^ 3 ]. The critical tempera- 
ture and the Pippard coherent length are given by 



T c = —A w 0.57A; 

7T 
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(22) 



This completes the formalism of formation of Cooper 
pairs with Majorana neutrinos and BCS superconduc- 
tivity. This may have many applications. 

We shall now discuss how the Majorana neutrino su- 
perfluidity, can solve the stability problem 0, [B[ in neu- 
trino dark energy models 0,0,0 ■ We demonstrate this in 
a specific two-generation neutrino dark energy model 0] • 
The standard model is extended with two right-handed 
neutrinos Ni, (i — 1, 2) and two scalars <&i, (i = 1, 2) with 
a global U(l)i x U(l)2 symmetry, so that these fields in- 
teract as 



C M = -aiJViJVJ^i + a 2 N 2 N c 2 <5> 2 . (23) 

When the fields <F; acquire vacuum expectation values 
(vev) ($,) = fi, we can express them as 



J_ 

2VT 



(24) 



where 4>i are the massless Nambu-Goldstone bosons and 
we assumed same decay constant fi — f for both the 
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fields. Writing = Mjexp[2i<^//], we get the masses 
of the right-handed neutrinos Ni to be Mj. 

In this model, the neutrino Dirac mass terms rriijViNj 
do not respect the global symmetry, and hence, one com- 
bination of the global U(l) symmetries is broken ex- 
plicitly. As a result, one of the two massless Nambu- 
Goldstone bosons picks up a small mass, making it a 
pseudo Nambu-Goldstone boson (pNGB). This pNGB 
(denoted by <j>) can then become the acceleron field that 
explains the dark energy 

After integrating out the heavy right-handed neutri- 
nos, we can write down the mass matrix of the light 
physical neutrinos vE = ( v\ v 2 ) , as given in equation 
SI where = —im/2Mf is not diagonal, leading to a 
long-range attractive force, and m p = m 2 /M. In these 
models of MaVaNs, naturalness restricts the mass scales 
of the model and the right-handed neutrino mass scale is 
supposed to be as low as eV. 

To estimate the parameter x we need to know C which 
depends on the acceleron mass C ~ ahr- If °ne takes 
rritj, ~ 6Mpc [5(, rn v ~ leV and n v ~ 56/cm 3 one finds 
x ~ 1CP 56 . Thus practically the exponential factors in 
Eqs. (21-23) are unity. So, the relevent scales of this equa- 
tions are determined by the cut-off A and neutrino mass 
m v . There remains a great deal of uncertainty over the 
range of these parameters. For example if one can take A 
as the scale when the tracking p v ~ Pd_b becomes valid 
[see Fardon et.al in Ref.[2j], one can take A as the de- 
coupling temperature 1 MeV. One can also take A to be 
very close to the scale when the neutrinos become non- 
relativistic i.e. few times m v . If we take the neutrino 
mass around 1 eV, £ has range between 0.36 — 10 4 cms. 
The values of £ in this entire range still can be smaller 
than the neutrino lumps [l8T ] . Finally we comment on the 
so called instability in MaVaN models Q in this changed 
scenario with the condensates. The instability is known 
to arise when the coupling between the scalar field and 
the neutrinos is stronger than the gravitational force. 
The coupling can be described by a source term of the 
type 0{(j)){p v — ip v ) [5]. In the the relativistic regime 
p v ~ 3p Ul the coupling is highly suppressed. However, it 
can become strong in the non-relativistic limit p v ~ 0. 
Size of the Cooper pairs is determined by the interaction 
strength between the scalar field and the neutrinos. <f> s 
dynamics is given by the following Lagrangian density 



C = d <fitd <f> s - di4>ldi(j) s - V{(j> s ), 



(25) 



where, m ~ 2to„ represents mass of the condensate. The 
condensate potential 



V(ct> s )=m 2 \<j> s \ 2 + g\<f> s \ 4 



(26) 



There can be interaction between the condensate </> s and 
the scalar field <p given by V\ n t = 9\4>4- 's- We can write 
down the perturbation equation for the coupled system 



4> and <fis, following Bjaelde, et. al. in Ref.[5], for the 
minimum of the effective potential tracked by the fields, 



! + 2H5$ + [k 2 + a 2 V^'] 
2H64> S + [k 2 + a 2 V^' s ] S 



= giHs (27) 

= giH- (28) 



It should be noted that g\ quantifies the coupling be- 
tween the condensate field <p s and the scalar field <f> and 
it can be much smaller than the coupling between the 
neutrinos and the scalar field. The condensates already 
formed due to the attractive between the Majorana neu- 
trino mediated by the <fi. 

Eqs. (|27ll25)) contain VI' and Vi' terms which are non- 
linear functions of <p and (j) s respectively. These equations 
are linearized by making the assumption that they can 
be written as (f> ~ (f>o + 5<fi and <p s w 4> s q + 8<j> s , where, 
the quantities with suffix represent background quan- 
tities and they can vary with time much slowly than the 
perturbation. From this perturbation in the scalar field 



gi8(j) s 



2iH + k 2 + a 2 V" 



(29) 



We can compute VI in Eq.(24) using two forms of po- 
tential Vj, frequently used in the MaVaN models. First 
we consider Coleman- Weinberg type of potential [n| 



V4, = V Q log(l + K(j>) 



(30) 



where parameters Vq and k can be selected to yield to- 
tal dark energy contribution Qde — 0.7. Using the 
above linearization one can find V'l = , 2 ^ y% . Secondly 
we consider an inverse power-law model of the potential 
given by 



Vi> 



AI 



n+4 



(31) 



where the parameter M can be fixed by the requirements 
Qde — 0.7 and 3> H. Thus we write as for n < 1, 



VI = 



n(n - l)(n-2)M 



n+4 



(32) 



Instability is known to occur in a MaVaN scenario for 
H < k/a < uifp [||. The following dispersion relation can 
be obtained from Eqs. (|27H28p : 

ui 4 - 4iHu; 3 - [AH 2 + B] lu 2 + 2iHBuj + B X B 2 - g\ = 

(33) 
and 



where, B = 2k 



B 2 = k 2 



+ « 2 (Vf + V£), B, = k 2 + a 2 V{' 



a V'!. Eq.(30) is quartic in uj and it can 



be solved exactly by analytical means. However it is in- 
structive to solve it by approximate methods. The terms 
that are linear in the expansion rates H will contribute 
to the damping of the modes. Moreover the instability, in 
the MaVaN scenario without the condensates, the insta- 
bility was found to occur in the regime H < k/a < m^. 
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Thus H defines the smallest wave- vector of the perturba- 
tions. In what follows we ignore the terms involving H in 
Eq. (f3"3")) to get following biquadratic dispersion relation, 



lo 4 - Blu 2 + B\B 2 - gf « 
Solution of this can be written as 



= B ± JB 2 -A{B 1 B 2 



9l) 



(34) 



(35) 



None of the roots of Eq.(32) is imaginary if B > and 
BiB 2 — g\ > 0. These are the stability criteria and their 
validity may depend on the form of the scalar field po- 
tential as B and Bi both involve the term VI'. If B > 
and BiB 2 ~ g\ < then the two roots of Eq.(39) are 
purely imaginary and one of it could give the instabil- 
ity. It is easy to verify that for the Coleman- Weinberg 
potential given by Eq.(33) V£ is positive for the parme- 
ters values V » 8.6 x 10" 13 e^ 4 , «» lx lO^M^ 1 and 
(f)Q w 10~ 6 Mpi taken from Ref. For the condensate 
potential with g > 0, one can have all the stability con- 
ditions satisfied for a sufficiently small coupling strength 
g% . The stability condition satisfied for the entire regime 
H < k/a < 77i0. For the power-law kind of potential 
given by Eq.(34) one can take M sa O.OlleV, n = 0.01 
and O « 0.001Mp Z @ ,one can have B\ > 0. For this 
case the stability conditions are again satisfied for a suf- 
ficiently small values of the coupling g\ and g > 0. 

The Majorana neutrino superfluidity we discussed may 
have some interesting consequences. The effects of the 
attractive long force, required for the formation of con- 
densates, have been discussed for the neutrino oscillation 
experiments and also in cosmology [3] ■ The acceleron 
potential can also change some of the features of the neu- 
trino oscillations 2l|, which will be further modified in 
the presence of the condensates. 

We would like to note here that the stability analysis 
provided in this paper is for the scalar field 4> s describing 
neutrino condensates. This is a valid description when 
the Bose-Einstein condensation (BEC) is formed. How- 
ever when the coherent length is much smaller than the 
neutrino interparticle spacing, the phenomenon of super- 
fluidity of Majorana neutrino still occur in BCS as shown 
in Ref.([l[). Furthermore, the accelerated expansion can 
make the ratio of the helicity flip rate to the Hubble ex- 
pansion rate even more smaller than the one we have at 
present. Thus the condition for defining chemical poten- 
tial in accelerated universe will remain satisfied in future. 
However the chemical potential thus defined is a 'dynam- 
ical quantity'. But the question of the small variation in 
the chemical potential with time may not be studied by 
the formalism presented here. 

In summary we have discussed two issues for MaVaN 
models namely the superfluidity of Majorana neutrinos 
and the stability of MaVaN dyanmics. We proposed a 
formalism to have condensates with Majorana neutrinos. 



Our formalism shows that the neutrino superfluidity nat- 
urally arises in MaVaN scenario and is a generic feature 
of interacting Fermi particles having attractive potential 
between them at a low temperature (22| . In addition, we 
have also shown that for the case when the condensate 
dynamics become important and > H the dynamics 
of the condensate can be stable for a variety of the dark 
energy potentials. 
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